MCMC analysis of WMAP3 and SDSS data points to broken symmetry inflaton 
potentials and provides a lower bound on the tensor to scalar ratio. 
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We perform a MCMC (Monte Carlo Markov Chains) analysis of the available CMB and LSS 
data (including the three years WMAP data) with single field slow-roll new inflation and chaotic 
inflation models. We do this within our approach to inflation as an effective field theory in the 
Ginsburg-Landau spirit with fourth degree trinomial potentials in the inflaton field (f>. We derive 
explicit formulae and study in detail the spectral index n s of the adiabatic fluctuations, the ratio r 
(N of tensor to scalar fluctuations and the running index dn„/dlnk. We use these analytic formulas as 

(— ( ■ hard constraints on n 3 and r in the MCMC analysis. Our analysis differs in this crucial aspect from 

previous MCMC studies in the literature involving the WMAP3 data. Our results are as follow: 
(i) The data strongly indicate the breaking (whether spontaneous or explicit) of the <f> — > — <j> 
• symmetry of the inflaton potentials both for new and for chaotic inflation, (ii) Trinomial new 

inflation naturally satisfies this requirement and provides an excellent fit to the data, (iii) Trinomial 
chaotic inflation produces the best fit in a very narrow corner of the parameter space, (iv) The 
^ • chaotic symmetric trinomial potential is almost certainly ruled out (at 95%CL). In trinomial 

' chaotic inflation the MCMC runs go towards a potential in the boundary of the parameter space 

and which ressembles a spontaneously symmetry broken potential of new inflation, (v) The above 
results and further physical analysis here lead us to conclude that new inflation gives the best 
description of the data, (vi) We find a lower bound for r within trinomial new inflation potentials: 
r > 0.016 (95% CL) and r > 0.049 (68% CL). (vii) The preferred new inflation trinomial potential is 
a double well, even function of the field with a moderate quartic coupling yielding as most probable 
values: n s ~ 0.958, r ~ 0.055. This value for r is within reach of forthcoming CMB observations. 
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I. INTRODUCTION AND RESULTS 

Inflation was introduced to solve several outstanding problems of the standard Big Bang model [l[ and has now 
become an important part of the standard cosmology. At the same time, it provides a natural mechanism for the 
generation of scalar density fluctuations that seed large scale structure, thus explaining the origin of the tempera- 
ture anisotropics in the cosmic microwave background (CMB), as well as that of tensor perturbations (primordial 
gravitational waves) @, 0, 0] . 

A distinct aspect of inflationary perturbations is that these are generated by quantum fluctuations of the scalar 
field(s) that drive inflation. After their wavelength becomes larger than the Hubble radius, these fluctuations are 
amplified and grow, becoming classical and decoupling from causal microphysical processes. Upon re-entering the 
horizon, during the matter era, these classical perturbations seed the inhomogeneities which generate structure upon 
gravitational collapse 0, A great diversity of inflationary models predict fairly generic features: a gaussian, nearly 
scale invariant spectrum of (mostly) adiabatic scalar and tensor primordial fluctuations, making the inflationary 
paradigm fairly robust. The gaussian, adiabatic and nearly scale invariant spectrum of primordial fluctuations provide 
an excellent fit to the highly precise wealth of data provided by the Wilkinson Microwave Anisotropy Probe (WMAP) [f| 
Perhaps the most striking validation of inflation as a mechanism for generating superhorizon ('acausal') fluctuations 
is the anticorrelation peak in the temperature-polarization (TE) angular power spectrum at I ~ 150 corresponding 
to superhorizon scales [5[ . The confirmation of many of the robust predictions of inflation by current high precision 
observations places inflationary cosmology on solid grounds. 

Amongst the wide variety of inflationary scenarios, single field slow roll models provide an appealing, simple and 
fairly generic description of inflation. Its simplest implementation is based on a scalar field (the inflaton) whose 
homogeneous expectation value drives the dynamics of the scale factor, plus small quantum fluctuations. The inflaton 
potential is fairly flat during inflation. This flatness not only leads to a slowly varying Hubble parameter, hence 
ensuring a sufficient number of e-folds, but also provides an explanation for the gaussianity of the fluctuations as well 
as for the (almost) scale invariance of their power spectrum. A flat potential precludes large non-linearities in the 
dynamics of the fluctuations of the scalar field. 

The current WMAP data seem to validate the simpler one- field slow roll scenario 0, Q. Furthermore, because the 
potential is flat the scalar field is almost massless, and modes cross the horizon with an amplitude proportional to 
the Hubble parameter. This fact combined with a slowly varying Hubble parameter yields an almost scale invariant 
primordial power spectrum. The slow-roll approximation has been recently cast as a 1 /N expansion Q , where N ~ 50 
is the number of efolds before the end of inflation when modes of cosmological relevance today first crossed the Hubble 
radius. 

The observational progress permits to start to discriminate among different inflationary models, placing stringent 
constraints on them. The upper bound on the ratio r of tensor to scalar fluctuations obtained by WMAP 0, Q 
necessarily implies the presence of a mass term in the single field inflaton potential and therefore rules out the 
massless monomial (j) 4 potential @, 0). Hence, as minimal single field model, one should consider a sufficiently general 
quartic polynomial, that is the trinomial potential. 

Besides its simplicity, the trinomial potential is a physically well motivated pot ential for inflation in the context of 
the Ginsburg-Landau approach to effective field theories (see for example ref.llOj]). This potential is rich enough to 
describe the physics of inflation and accurately reproduce the WMAP data 

The slow-roll expansion plus the WMAP data constraints the inflaton potential to have the form Q 

V{4>) = N M 4 w(x) ! (1.1) 
where 4> is the inflaton field, x is a dimcnsionless, slowly varying field 

(1.2) 



N Mpi 
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w(x) ~ O(l) and M is the energy scale of inflation which is determined by the amplitude of the scalar adiabatic 
fluctuations @ to be 

M ~ 0.003f9 Mpi = 0.77 x 10 16 GeV . (1.3) 

Following the spirit of the Ginsburg-Landau theory of phase transitions, the simplest choice is a quartic trinomial for 
the inflaton potential @, H, Q : 

»W = ».±ix 2 + f\/fx» + |x«, d.4) 

where the coefficients wo, h and y are dimensionless and of order one and the signs ± correspond to large field 
and small field inflation, respectively (namely, chaotic inflation and new inflation, respectively), h describes how 
asymmetric is the potential, y is the dimensionless quartic coupling. 
Inserting eq. (|1.4[) in cq. ljl.ip yields, 

^) = F ±^V + ^V + ^ 4 , (1.5) 
where the mass m 2 and the couplings g and A are given by the following see-saw-like relations, 

2 / > ,r \ 4 



Notice that y ~ 0(1) ~ h guarantee that g ~ 0(1O~ 6 ) and A ~ 0(1O~ 12 ) without any fine tuning as stressed in ref. 
Q • That is, the smallness of the couplings here is not a consequence of fine tunning but follow directly from the form 
of the inflaton potential eq. (|l.l[) and the amplitude of the scalar fluctuations that fixes M [7J]. 

During inflation the inflaton field x slowly runs from its initial value till its final value Xend at the absolute minima 
of the potential itf(x). We have Xend = for chaotic inflation while Xend turns to be a function of the coupling y and 
the asymmetry h for new inflation. Inflation ends after a finite number of efolds provided 

w{Xend) = w'{Xend) = 

Enforcing this condition in the inflationary potential eq. (|1.4p determines the constant wq. We have wo = for chaotic 
inflation while wq turns to be a function of the coupling y and the asymmetry h for new inflation, x can vary in the 
interval (0, oo) for chaotic inflation while x is m th e interval (0, Xend) for new inflation. 

We derive explicit formulae and study in detail the spectral index n s of the adiabatic fluctuations, the ratio r of 
tensor to scalar fluctuations and the running index dn s /d\nk both for trinomial chaotic inflation and trinomial new 
inflation. 

The small coupling limit y — > of eqs. (jl.4[) - (|1.5[) corresponds to the quadratic monomial potential while the strong 
coupling limit y — > oo yields in chaotic inflation the massless quartic monomial potential. The extreme asymmetric 
limit \h\ — > oo yields a massive model without quadratic term. In such a limit the product \h\ M 2 must be kept fixed 
since it is determined by the amplitude of the scalar fluctuations. 

In this paper we perform Monte Carlo Markov Chains (MCMC) analysis of the commonly available CMB and Large 
Scale Structure (LSS) data. For CMB we considered the three years WMAP data, which provide the dominating 
contribution, and also small scale data (ACBAR, CBI2, BOOMERANG03). For LSS we considered SDSS (Sloan 
Digital Sky Survey) [ll|. We used the CosmoMC program [l2j within the effective field theory of inflation. We used 
a large collection of parallel chains with a total number of samples close to five million. For chaotic and new inflation 
we imposed as a hard constraint that the spectral index n s of the adiabatic fluctuations and the ratio r are given 
by the analytic formulas at order 1/N we derived for the trinomial inflaton potential. Our analysis differs in this 
crucial aspect from previous MCMC studies involving the WMAP3 data set [11]. As natural within inflation, we 
also included the inflationary consistency relation nx — —r/8 on the tensor spectral index. This constraint is in any 
case practically negligible. The details of the MCMC analysis are explained in Section V. 

Our approach is different to the inflationary flow equations [l4[ where the inflaton potential changes (that is, the 
model changes) as the flow goes on. We work with a given potential within the Ginsburg-Landau (GL) spirit, that 
is the trinomial potential. We investigate the physics of the chosen potential in the parameter space driven by the 
data through the Monte Carlo Markov Chains. In our work, n s and r arc computed analytically to order 1/N for the 
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trinomial potentials [eas. i|3.6| )- (|3.8|) and (|4.6|) - (|4.8p ]. Since N ~ 50, higher order corrections in 1/N are irrelevant and 
can be safely neglected. As shown in ref. [7j, the slow roll expansion is in fact a systematic expansion in 1/N. 

We allowed seven cosmological parameters to vary in our MCMC runs: the scalar spectral index n s , the tensor-scalar 
ratio r, the baryonic matter fraction u>b, the dark matter fraction ui c , the optical depth r, the ratio of the (approximate) 
sound horizon to the angular diameter distance 9 and the primordial superhorizon power in the curvature perturbation 
at 0.05 Mpc" 1 , A s . We allowed the same seven parameters to vary in the MCMC runs for chaotic and new inflation. 

In the case of new inflation, since the characteristic banana-shaped allowed region in the (n s ,r) plane [fig. d] is 
quite narrow and non-trivial, it is convenient to use the two independent variables z = | x 2 an d h in trinomial 
inflationary setup as MC parameters. That is, we used the analytic expressions we found at order 1/N, to express n s 
and r in terms of z and h. 

Concerning priors, we kept the same, standard ones, of the ACDM+r model for the first five parameters u) c , t, 8 
and A s ), while we considered all the possibilities for z and h. 

In the case of chaotic inflation we kept n s and r as MC parameters, imposing as hard priors that they lay in 
the region described by chaotic inflation [fig. 2]. This is technically convenient, since this region covers the major 
part of the probability support of n s and r in the ACDM+r and the parametrization eqs. (|3.6p - (|3.8p in terms of the 
parameters z and h becomes quite singular in the limit h — > — 1. This is indeed the limit which allows to cover the 
region of highest likelihood. The priors on the other parameters where the same of the ACDM+r model and of new 
inflation. 

We did not marginalize over the SZ amplitude, and we did not include non-linear effects in the evolution of the 
matter spectrum. The relative corrections are in any case not significant [5], especially in the present context. 

In all our MCMC runs we keep fixed the number of efolds N since horizon exit till the end of inflation. The reason is 
that the main physics that determines the value of N is not contained in the available data but involves the reheating 
era. Therefore, it is not reliable to fit N solely to the CMB and LSS data. The precise value of N is certainly near 
N = 50 [14, 15]. We take here the value N — 50 as a reference baseline value for numerical analysis, but from the 
explicit expressions obtained in the slow roll 1/N expansion we see that both n s — 1 and r scale as 1/N. Therefore, 
varying N produces a scale trasformation in the (n r — 1, r) plane, thus displacing the black and red curves in fig. [2] 
towards up and left or towards down and right. This produces however small quantitative changes in our bounds for 
r as well as in the most probable values for r and n s . MCMC simulations with variable N and imposing the trinomial 
new inflation potential yielded N ~ 50 as the most probable value. 

We plot in fig. [2] r vs. n s both for chaotic and new inflation for fixed values of the asymmetry h and the 
coupling y varying along the curves. We see that the regions of the trinomial new inflation and chaotic inflation are 
complementary in the (n s ,r) plane. The red curves are those of chaotic inflation, while black curves are for new 
inflation. The leftmost black line corresponds to new inflation with a symmetric potential h — 0. The rightmost black 
line describes the case of new inflation with an extreme asymmetric potential \h\ 3> 1. This last line is the border of 
the region on the right described by chaotic inflation. Although chaotic inflation covers a much wider area than new 
inflation, this wide area is only a small corner of the parameter space (field z = | \ 2 1 asymmetry h) as shown by fig. 

m 

For the trinomial new inflation model we find a lower bound on r: 

r > 0.016 (95% CL) , r > 0.049 (68% CL) (new inflation) . (1.7) 

while for n s we find: 

n s > 0.945 (95% CL) (new inflation) . 

The most probable values are (see fig. |6]), 

n s ~ 0.956 , r ~ 0.055 (new inflation) . (1.8) 

For trinomial new inflation there exists the theoretical upper limits: n s < 0.9615 . . . , r < 0.16 0. Thus, the most 
probable value of n s for trinomial new inflation eq. (|1.8p is very close to its theoretical limiting value, and that of r is 
not too far from it (see also fig. [6]). 

The probability distribution for the asymmetry parameter h is peaked at h = with 



\h\ < 4.92 with 95% CL new inflation 
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That is, we find that the most probable trinomial new inflation potential is symmetric and has a moderate nonlinearity 
with the quartic coupling y ~ 2 for h ~ 0. This is the following potential: 

The x ~~ * — X symmetry is here broken since the absolute minimum of the potential is at x 0- 

For trinomial chaotic inflation, the chaotic symmetric trinomial potential h = is almost certainly ruled out 
since h < —0.7 at 95 % confidence level (see fig. [SJ. 

We see the maximun probability for strong asymmetry h < —0.95 and significant nonlinearity 4.207 . . . < y < +oo 
That is, in chaotic inflation all three terms in the trinomial potential w(x) do contribute. 

We have not introduced the running of the spectral index dn s /dhik in our MCMC fits since the running must 
be very small of the order (D(N~ 2 ) ~ 0.001 in slow-roll and for generic potentials [7|. We find that adding the new 
parameter dn s /d\nk to the MCMC analysis make negligible changes on the fit of n s and r. All this suggests that the 
present data are not yet precise enough to allow a determination of dn s /dh\k. 

In summary, the data favour the breaking of the \ ^ ~X symmetry both for new and for chaotic inflation. 

Trinomial chaotic inflation produces its best fit to the data in a very narrow corner of the parameter space. 

The potential in chaotic inflation has a single minimum. But this minimum for the best fit potential happens to 
be in the boundary of the parameter space, precisely where the asymmetry of the potential is so large that arises 
an extra minimum to the potential. That is, the data are begging the potential to be a double well (that is, a new 
inflation potential) . The MCMC runs go towards a potential in the boundary of the parameter space and maximal 
symmetry breaking. This limiting potential exhibits an inflexion point and a significative nonlinearity. This strongly 
suggest that the true potential may be outside the parameter space of chaotic inflation. On the contrary, starting 
with a double well, that is new inflation, the data are perfectly well fitted with small or even zero asymmetry. 

In chaotic inflation, the MCMC minimization leads to a large nonlinearity (from the higher order cubic and cuartic 
terms), making chaotic inflation not stable under higher order terms in the potential. Chaotic inflation thus results 
contrary to the Ginsburg-Landau spirit in which higher order terms of the potential must be smaller and smaller 
giving irrelevant corrections. 

On the contrary, new inflation results perfectly well fitted by the data with very small or zero cubic non-linearity, 
and thus new inflation is stable under higher order terms of the potential, which is perfectly well within the Ginsburg- 
Landau sprit. 

Indeed, one can fit the WMAP+LSS data both with trinomial new inflation as well as with trinomial chaotic 
inflation. However, the trinomial potential that gives the best fit for chaotic inflation is theoretically undesirable 
for the reasons explained above. Therefore, trinomial new inflation is the best model that describes the data. This 
implies that we do have a lower bound on r as given by eq. (|1.7[) . Moreover, the best fit is obtained with the double 
well inflation potential eg. (|1 .9|) yielding the values in cq. (jl.8p for n s and r. 

This paper is organized as follows: in section II we discuss the slow-roll approach as the l/N expansion. In sections 
III and IV we present the trinomial chaotic inflation and the new trinomial inflation, respectively, displaying the 
analytic formulas for the spectral index n s of the adiabatic fluctuations, the ratio r of tensor to scalar fluctuations 
and the running index dn s /d\nk. In section V we explain our MCMC analysis, we present our MCMC results for 
new and chaotic trinomial inflation and we state our conclusions. 



II. THE INFLATON POTENTIAL AND THE l/N SLOW ROLL EXPANSION 

The description of cosmological inflation is based on an isotropic and homogeneous geometry, which assuming flat 
spatial sections is determined by the invariant distance 

ds 2 =dt 2 -a 2 (t) dx 2 . (2.1) 

The scale factor obeys the Friedman equation 
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where pit) is the total energy density and Mpi = 1/^/SnG = 2.4 x 10 18 GeV. 

In single field inflation the energy density is dominated by a homogeneous scalar condensate, the inflaton, whose 
dynamics is described by an effective Lagrangian 



C = a 3 it) 



^ (V0) 2 



2 a 2 it) 



Vi4>) 



(2.3) 



The inflaton potential V(d>) is a slowly varying function of <fi in order to permit a slow-roll solution for the inflaton field 
<Pit). We showed in ref.Q that combining the slow roll expansion with the WMAP data yields an inflaton potential 
of the form 



V(<j>) = N M 4 wix) 



where % is a dimensionless, slowly varying field 



X = 



N M 



(2.4) 



(2.5) 



pi 



wix) ~ 0(1) , iV ~ 50 is the number of efolds since the cosmologically relevant modes exited the horizon till the end 
of inflation and M is the energy scale of inflation 

The dynamics of the rescaled field x exhibits the slow time evolution in terms of the stretched dimensionless time 
variable, 



t M 1 



Mpi ViV 



n = 



H 



N m 



= 0(1) 



(2.6) 



The rescaled variables x an< l r change slowly with time. A large change in the field amplitude <p results in a small 
change in the x amplitude, a change in (f> ~ Mpi results in a x change ~ 1/yN. The form of the potential, eq. ()2.4j) 
and the rescaled dimensionless inflaton field eq. (|2.5D and time variable r make manifest the slow-roll expansion as a 
consistent systematic expansion in powers of 1/N |7fl. 

The inflaton mass around the minimum is given by a see-saw formula 



M 2 



M. 



2.45 x 10 13 GeV 



pi 



The Hubble parameter when the cosmologically relevant modes exit the horizon is given by 

H = VN mH ~ 1.0 x 10 14 GeV = 4.1 m , 

where we used that 7i ~ 1. As a result, m <C M and H <C Mpi. 

A Ginsburg-Landau realization of the inflationary potential that fits the amplitude of the CMB anisotropy remark- 
ably well, reveals that the Hubble parameter, the inflaton mass and non-linear couplings are see-saw-like, namely 
powers of the ratio M/Mpi multiplied by further powers of 1/N. Therefore, the smallness of the couplings is not 
a result of fine tuning but a natural consequence of the form of the potential and the validity of the effective field 
theory description and slow roll. The quantum expansion in loops is therefore a double expansion on iH/Mp{) 2 and 
1/N. Notice that graviton corrections are also at least of order iH/Mpi) 2 because the amplitude of tensor modes is 
of order H/Mpi . We showed that the form of the potential which fits the WMAP data and is consistent with slow 
roll eas. (|2.4[) - (l2.5[l implies the small values for the inflaton self-couplings 0]. 

The equations of motion in terms of the dimensionless rescaled field x an d the slow time variable r take the form, 



U\r) = 
N dr 2 



1 1 (d X 
2>y2N \dr 

^n^ + w'ix) 



wix) 



(2.7) 



The slow-roll approximation follows by neglecting the terms in eas. (|2.7p . Both wix) and 7Y(t) are of order N° for 
large N. Both equations make manifest the slow roll expansion as an expansion in 1/N. 
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The number of e-folds N[x] since the field x exits the horizon till the end of inflation (where \ takes the value Xend) 
can be computed in close form from eqs. (|2.7p in the slow-roll approximation (neglecting 1/N corrections) 



N[ X ] 
N 



w'(x) 



d X < 1 



(2.8) 



Inflation ends after a finite number of efolds provided 

w{Xend) = w'{Xend) = 

This condition will be enforced in the inflationary potentials. 

The amplitude of adiabatic scalar perturbations is expressed as 0, 0, M, 0] 

, A (5) |2 = _n^_ [jlY w3 m 

1 kadl 12 tt 2 \Mpi) W '\x) ' 
The spectral index n s , its running and the ratio of tensor to scalar fluctuations are expressed as 



(2.9) 



(2.10) 



n H — 1 = 

N 



w'{x) 



i 2 



w(x) 



2_ w"( X ) 
N w(x) 



dn s 
dink 



N 



2 w'( X )w'"(x) 6 K( X )] 4 + 8 [w'(x)?w"(x) 



N 2 w 2 (x) 
w'(x) 



N 2 w 4 (x) N 2 w 3 (x) 



(2.11) 



In eas. (|2.8p - (|2.1ip the field x is computed at horizon exiting. We choose N[x] = N = 50. 
Since, w(x) and w'{x) are °f order one, we find from eq. (|2.10[) 



/ M 



pi 



^|Af) |^1.02xl0- 5 

I A: ad I 



(2.12) 



where we used N ~ 50 and the WMAP value for |A^ d | = (4.67 ± 0.27) x 10~ 5 5]. This fixes the scale of inflation 
to be 

M ~ 3.19 x 10 -3 A/pl =i 0.77 x 10 16 GeV . 

This value pinpoints the scale of the potential during inflation to be at the GUT scale suggesting a deep connection 
between inflation and the physics at the GUT scale in cosmological space-time. 

We see that \n s — 1| as well as the ratio r turn out to be of order 1/A^. This nearly scale invariance is a natural 
property of inflation which is described by a quasi-de Sitter space-time geometry. This can be understood intuitively 
as follows: the geometry of the universe is scale invariant during de Sitter stage since the metric takes in conformal 
time the form 



ds 2 



1 



[(d V f - (dx) 



Therefore, the primordial power generated is scale invariant except for the fact that inflation is not eternal and lasts 
for TV efolds. Hence, the primordial spectrum is scale invariant up to 1 /N corrections. The Harrison-Zeldovich values 
n s = 1, r = and dn s /d\n.k = correspond to a critical point as discussed in ref.[7[. This gaussian fixed point is 
not the inflationary model that reproduces the data but the inflation model hovers around it in the renormalization 
group sense with an almost scale invariant spectrum of scalar fluctuations during the slow roll stage. 

We analyze in the subsequent sections chaotic inflation and new inflation in its simple physical realizations within 
the Ginzburg-Landau approach (the trinomial potential) [7J, |g, l9( . We perform with these models Monte Carlo Markov 
chains analysis of the three years WMAP data and LSS data. 
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(y/8) 1/2 x 

FIG. 1: Trinomial Chaotic Inflation. We plot here the chaotic inflation trinomial potential [eg. (|3.2p with positive quadratic 
term] | w(x) vs - the field variable y/z = \/y/8 x f° r different values of the asymmetry parameter h, namely, h = 0, —0.5, —0.8 
and —1. Notice the inflection point at y/z — 1 when h = — 1. 



III. TRINOMIAL CHAOTIC INFLATION: SPECTRAL INDEX, AMPLITUDE RATIO, RUNNING 

INDEX AND LIMITING CASES 

We consider now the trinomial potential with unbroken symmetry investigated in ref. Q : 

V((p) = — (p + — ( P +4^ . (3-1) 
here m 2 > and g and A are dimensionless couplings. 

The corresponding dimensionless potential w(x) eas. (|l.l|l - (|1.2l ) takes the form 

1 -,a , h [V „3 , 2/ „4 



-W = 2 X + 3Vf X + ^ X ' (3 - 2) 

where the quartic coupling y is dimensionless as well as the asymmetry parameter h. The couplings in eq. (|3.ip and 
eq. (|3.2|) are related by eq. (|1.6p . 

Chaotic inflation is obtained by choosing the initial field x m the interval (0, +oo). The inflaton x slowly rolls down 
the slope of the potential from its initial value till the absolute minimum of the potential at the origin. 

Computing the number of efolds from eq. (|2.8p . we find the field x a t horizon crossing related to the couplings y 
and h. Without loss of generality, we choose h < and shall work with positive fields x- 

The potential eq. (|3.2p has extrema at x = and %± given by, 



X± = y ~ [~h ± i A] , A= y/l - h 2 . (3.3) 

That is, for \h\ < 1, w(x) has only one real extremum at x = while for \h\ > 1, w(x) has three real extrema. There 
is always a minimum at x = since w/'(0) = 1. At the non-zero extrema we have 

w"(x±) = -2 A (A±ih) . 
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We have for \h\ > 1, 

and w"(x±) = 2\Jh 2 -\ (V 'h 2 -l^fc) . 

Hence, for any /i < — 1, we have w"(x+) > while «/'(%_) < 0. Notice that x± > for h < —1. 

Therefore, we have chaotic inflation for positive field x m the regime \h\ < 1 using the inflaton potential ea. (|3.2|) . 

We can also have chaotic inflation with the potential eq. (|3.2[) for negative field if h > 3/V8, but for \h\ > 3/v8 the 
absolute minimum is no more at x — but at x+- Since w(x+) < for \h\ > 3/^/8 we have to subtract in this case 
the value w(x+) from w(x) in order to enforce ea. (|2.9|) . 

We consider in subsections IIII A l IIII Bl and IIII CI the regime — 1 < h < 0, x > 0. The case h < —1 is analyzed in 
subsection IIII D[ 



X± 



-h±\Jh 2 - 1 



A. The small asymmetry regime: — 1 < h < 0. 



In chaotic inflation the inflaton field slowly rolls down the slope of the potential from its initial value till the absolute 
minimum of the potential at x = 0. 



It is convenient to define the field variable z by: 



_ v 2 

Z = o x 



(3.4) 



In terms of z the chaotic trinomial potential takes the form 

, . 4 / 4, _ 1 
W (X) = - z I 1 + - ft V z + 2 z 

When < 1 we are in the quadratic regime where w(x) is approximated by the x 2 term. For z > 1 we go to the 
non-linear regime in z and all three terms in w(x) are of the same order of magnitude. 

In fig. [1] we plot | w(x) as a function of yfz for several values of h > — 1. We see that the potential becomes flatter 
as h decreases. For h = — 1, both w'(x) an d w'(x) vanish at y/z = 1. The case h = —1 is singular since the inflaton 
gets stuck an infinite amount of time at the point \fz = 1. 

By inserting eq. (|3.2[) for w(x) into ea. (|2.8[) for iV[x] and setting N[x] — N we obtain the field x or equivalently z , 
at horizon exit, in terms of the coupling y and the asymmetry parameter h, 



y 



z + -h^+[l--h 2 \ log (l + 2 h sfz + z) 



Ah 
3A 



2h' 



arctan 



A 



arctan 



(3.5) 



This defines the field z as a monotonically increasing function of the coupling y for < y, z < +oo. Recall that x 
and z corresponds to the time of horizon exit. 

We obtain from eqs. (|2~T0l) . (|2TTT|) and |3~2|) the spectral index, its running, the ratio r and the amplitude of adiabatic 
perturbations, 



n s = 1 — 



2iVz 



(l + 2/i yfl+z) 



l + 4hyfz + 3z 



(3.6) 
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2z 2 TV 2 



(l + 2/l y/z~+ Z )(h y/z+\z) (l + 2h^+zf 

■ O 



(1 + 2 h y/z + z) (1 + 4 fe + 3 z) 
(l + i^Vi+iz) 3 



(3.7) 



r = 



4y (1 + 2^V^+^) 2 



(3.8) 
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( S) 2 = 2TV 2 f M\ 4 z 2 (l + lh^+^zf 



(3.9) 



In chaotic inflation, the limit z — » + implies y — > + (shallow limit), we have in this limit: 

y 2z + 0(z?) 



= 1 77 

TV 



8 (s) 2 ,-0+ A^_ f M \ 

TV ' I * arfl 



6tt 2 Im p J 



(3.10) 



The results in the y — > + limit are independent of the asymmetry ft and coincide with those for the purely quadratic 
I 

2 



monomial potential i x 2 . 



In the limit z — > +oo which implies y — > +oo (steep limit), we have for fixed ft > — 1, 



z + -/i^+ l--/i 2 ] logz + C(l) 



y— *+oo 

1G 

TV 

.(S) 1 2 y^_+oo 
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4 ft 
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■oil 
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4 ft 

3 7^ 

N 2 ( M 



1 



12 7T 2 
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For ft = 0, n s and r in the limit y — > +00 coincide with those of the purely quartic monomial potential i x 4 ' 



n s = 1 

N 



16 

TV 



(3.11) 



(3.12) 



B. The flat potential limit h 



Namely, w'(xo) — w"(xo) = while 



We consider chaotic inflation in the regime — 1 < ft < 0, x > 0. 

At ft = —1 the potential eq. (|3.2|) exhibits an inflexion point at xo = 

w(xo) = 2/(3 y) > 0. That is, this happens at zo — | Xo — 1- 

Therefore, for ft > — 1 but very close to ft = — 1, the field evolution strongly slows down near the point x — Xo- This 
strong slow down shows up in the calculation of observables when the field x a t horizon crossing is x > XOi namely, 
for z > Zo = 1. For x < Xo, that is z < 1, the slow down of the field evolution will only appear if z ~ 1. Therefore, 
the limit ft = — 1 is singular since the inflaton field gets trapped at the point z = 1. 



Let us derive the ft — * — 1 + limit of y, n s and r from eqs. (|3.5[ )- (|3.9[) in the regimes z < 1 and z > 1, respectively. 

When ft — > — 1 + we see from eq. (l3.3[) that A — > and the arguments of the two arctan in eq. (|3.5[) diverge. Hence, 
the arctan tends to +7r/2 or —tt/2. Depending on whether z < 1 or z > 1, the 7r/2's terms cancel out or add, 
respectively. The special case z = 1 is investigated in the next subsection IIII CI 



In the case z < 1 we get: 
y = z 

n s = 1 - 



y 



2N z 



(l-Vi)(l-3Vi) 



/z+fz 



z < 1 , 

/!->-!+ 



z < 1 



(3.13) 
(3.14) 
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In particular, in the regime z — > 1~ 



(1 



y 2 (i - 



(3.15) 
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we find, 
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N 

2 1 

3 1 -VI 
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> o , /i 
AT M 2 



8tt Ml,; 
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(3.16) 



That is, in the limit h 



-1+, 



1 the ratio r becomes very small while the spectral index takes the value 

Xo = \ - because w'(xo) = and r is proportional 



n s = 0.92. The ratio r tends to zero in the regime z — > 1 , x 

to w' 2 (x) according to cq. (|2.11|l . 

The z — > 1 _ regime for /i — > — 1 + is a strong coupling limit since y — > +oo as shown by eq. (|3.16p . In addition, 
eq. (I3.16|) shows that for large y one must keep the product y M 2 fixed since it is determined by the amplitude of the 
adiabatic perturbations. We see from eq. (|3.16p that M = ^Jy M becomes the energy scale of inflation in the y — * oo 

limit: from cq. (|1.3|) . M ~ 10 16 GeV according to the observed value of |Aj^J/iV displayed in eq. ()2.12j) . while M and 
m vanish in the y —>■ oo limit 
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In fig. [2] we display r vs. n s for fixed values of the asymmetry parameter h and the coupling y varying along the 
curves. The red curves correspond to chaotic inflation. In the bottom of fig. H]we can see the curve for h — —0.999. 
The limiting curve for h — —1 (not drawn) will reach the point n s = 0.92 and the bottom line r = as described by 
ea. rf3~T6l) . 

We see from fig. [2] that the regions of the trinomial new inflation and chaotic inflation are complementary in the 
(n s , r) plane. New inflation describes the region of the (n s ,r) plane between the two black lines while chaotic inflation 
describes the whole plane to the right of the rightmost black line. 

In the case z > 1, we get from eas. (|3.5[) - (|3.9p . 

1]°) - +oo 



h + 1 
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z > 1 



z VhTT 
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\z) 2 



0{[h+l} ) 



<z+±z 



z > 1 . 



•0([ft + l] ) 



z > 1 , 



(3.17) 



In this strong coupling regime the indices become very large and hence in contradiction with the data. In addition, 
the slow roll expansion cannot be trusted when the coefficients of 1/N become large compared with unit. In conclusion, 
for h — > — 1 + the case of large field z > 1 is excluded by the data. 



C. The singular limit z = 1 and then h — » — 1 + yields the Harrison-Zeldovich spectrum 



We study in this section the case z — 1 in trinomial chaotic inflation. We obtain from eq. 



y^l + ih+(l-ihA logP(l + fc)]-!£(jj-2* a 



, h + l\ ( h 

arctan — - — — arctan — 
A J \A 



For z = 1, we find from eas.([3Tp ]l -(|3l ll and (|3~TS)l in the limit h -> -1+, 

1 



y 



=1. h-*- 1+ 7T / 1 

3 V 2{h + 1) ~ 3 



log(/i + l) + 0([/i+l] ) 



(3.18) 
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FIG. 2: Trinomial Inflation. We plot r vs. n s for fixed values of the asymmetry parameter h and the field z varying along the 
curves. The red curves are those of chaotic inflation with h < (only the short magenta curve has positive h), while the black 
curves are for new inflation. The color-filled areas correspond to 12%, 27%, 45%, 68% and 95% confidence levels according 
to the WMAP3 and Sloan data. The color of the areas goes from the darker to the lighter for increasing CL. New inflation 
only covers a narrow area between the black lines while chaotic inflation covers a much wider area but, as shown by fig. [5] 
this wide area is only a small corner of the field z - asymmetry h plane. Since new inflation covers the banana-shaped region 
between the black curves, we see from this figure that the most probable values of r are definitely non-zero within trinomial 
new inflation. Precise lower bounds for r are derived from MCMC in eq. (15.2[) . 



" 1 ' h = 1+ l + ^V2(ftTT) [l + O (/» + !)] , 



= 1, 16 7T 
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y/2(h + l)i[l + 0(h + l)] 



AS) ,2 z=l,h- 
*k ad\ 
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7r Mpi {h + !)■'■ 



(3.19) 



Therefore, we reach the Harrison-Zeldovich spectrum n s = 1, r = as a limiting value. This is a strong coupling 
regime y — > oo where in addition we must keep the ratio MJ (h + 1)' fixed for h — * — 1 + since it is determined by the 
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amplitude of the adiabatic perturbations. That is, we must keep 

M 



M 



- fixed as h — > — 1 + , 



(h + iy. 

while M as well as m go to zero. Actually, the whole potential V(<^>) eq. (|3.ip vanishes in this limit since, 



M 2 
W 



h— 1+ M 2 



pi 



M 



pi 



g = h 



nr ( M \ .=1,^-1+ nr / m y / i+j 



o, 



A = 



v ( M 



8 AT V Mpi 



=l, h-*-!" 1 



M 



24 N \Mpi 



1 + h 



(3.20) 



The inflaton field is therefore a massless free field at z — 1 and /i — > — 1 + . This explains why the corresponding 
spectrum is the scale invariant Harrison-Zeldovich one. This is clearly a singular limit since one cannot obtain any 
inflation from an identically zero potential. Namely, taking the flat limit h — > — 1 + in the spectrum computed for 
2 = 1, h > — 1 with a fixed number of efolds N, yields the scale invariant Harrison-Zeldovich spectrum. 

Notice that here we keep the number of efolds N fixed which makes the potential to vanish since otherwise the field 
will be stuck at the point 2 = 1 when h = — 1 leading to eternal inflation. 

In summary, this shows theoretically that the Harrison-Zeldovich spectrum is highly unplausible and unrealistic 
since it appears only in the singular limit z — 1, h — * — 1 + where the inflaton potential identically vanishes. 

Recall that one can also get a Harrison-Zeldovich spectrum letting formally the number of efolds TV to infinity in 
eqs. (f2~TT|) . 



D. The high asymmetry h < — 1 regime. 



In order to fulfill the finite number of efolds condition eq. (J2T9J) for h < —1 we have to add a constant piece to the 
chaotic inflationary potential eq. (|3.2p . 

We therefore consider as inflaton potential in the h < — 1 regime, 

The absolute minimum of this potential is at 

\- = x /^[-h + D] , D= yjh 2 1 , h<-l. (3.22) 



and we have 



G(h)= -h 4 -Ah 2 + l + -\h\D 3 



w 



u"( x+ ) = 2Vh 2 -l (^h 2 - 1 + > . (3.23) 

That is, the inflaton mass squared in units of m 2 takes the value 

2 D (D + \h\) . 

Here, the inflaton field rolls down the slope of the potential from its initial value larger than x+ till the absolute 
minimum of the potential at x — X+- 

By inserting eq. (l3.21|l for w(x) into eq. (|2.8p for N[x] and setting N[x] = N we obtain the field x or equivalently 
z = | x 2 j & t horizon exit, in terms of the coupling y and the asymmetry parameter h: 

y = z+^h^z + l + \h{D-h) + 2G{h) log ^l J - + ^h(h 2 -l)(D- h) ^g( ^ + 2 K D +D \ .(3.24) 
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This defines the field z as a monotonically increasing function of the coupling y for 

< y < +00 , z + = (D - h) 2 < z < +00 . 

Recall that \ and z corresponds to the time of horizon exit. 

We obtain from eas. (|2.10[) . (|2.1ip and (|3.21[) the spectral index, the ratio r and the amplitude of adiabatic pertur- 
bations, 



1 y 1 

n K = 1 



N (^/z'+h-D) 2 



Q z {^fz~ + h + D) 2 l + 4/iVi+3z 



[z + 2(D + !l)^-lh(D~h)~l] 2 z + 2 (D + §) ^ - § h (D - h) - 1 



16 V z {yfz + h + D) 2 

N (V^+h-D) 2 \ lZ + 2{D+ h ) ^_^ h{D _ h) _ l Y ' 

(S ) |2 _ JV^_ / M \ 4 [G{h) + 2z + l hz 3 / 2 + z 2 ] j^+h-D) 2 [z + 2(D + !p ^~lh{D-h)~l] 2 
12 7T 2 \Mpi) y 1 z (y/z + h + D) 2 



When yfz — > y'zZjT, y vanishes quadratically as, 

In this limit the spectral index, the ratio r and the amplitude of adiabatic perturbations become, 

y-o 2 8 (S ) pv-* N 2 f M 

s " TV ' N ' 1 feadl ~ 6vr 2 \M Pl 

These results are independent of the asymmetry h and coincide with those for the purely quadratic monomial 
potential \ \ 2 . 
We see here that 

8 16 3 2 

-<r<— , 1 - — < n s < 1 - — for < y < 00 . (3.26) 

Namely, the regime ft, < —1 of chaotic inflation covers values of r larger than the weak coupling limiting value r = ■£> 
eq. (|3.10p and smaller than the — > cxj pure quartic potential value r = 4f eq. (|3.12p . 

IV. TRINOMIAL NEW INFLATION: SPECTRAL INDEX, AMPLITUDE RATIO, RUNNING INDEX 

AND LIMITING CASES 

We consider here new inflation described by the trinomial potential with broken symmetry investigated in ref. [1, Q 

, , m 2 m a , A , , 

v(4>) = v - — <j> 2 + ^<t? + -<p\ (4.1) 

where m 2 > and g and A are dimensionless couplings. The corresponding dimensionless potential w(x) takes the 
form 

1 . 2 , h fy 3 y 4| 2 



-(x) = - 2 X< + ~ 3 ^ 2 X* + f 2 X* + ~F(h), (4.2) 

where the quartic coupling y is dimensionless as well as the asymmetry parameter h. The couplings in eq. (|4.ip and 
eq. (|4.2|) are related by, 

4 



and the constant wo [see eq. (|1.4|) ] is related to Vq in eq. ()4.ip by 

Wo = I F[h) = Ji^ 
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The constant F(h) ensures that w(x+) = w '{x+) — at the absolute minimum \ = \+ = y §(A + 1^1) 01 the 
potential w(x) according to eq. (|2.9p . Thus, inflation does not run eternally. F(h) is given by 

F(h) = | / l 4 + 4/ l 2 + 1 + || /i | A 3 ) A= \/h 2 + 1 . 

The parameter h reflects how asymmetric is the potential. Notice that w(x) is invariant under the changes % ~ > 
— x, h — > —h. Hence, we can restrict ourselves to a given sign for h. Without loss of generality, we choose h < and 
shall work with positive fields x as m sec. MIL 
We have near the absolute minimum x = X+i 

w( X ) X = X+ A(A + \h\) (x X+f + O (Vvlx - X+} 3 ) , (4.4) 
That is, the inflaton mass squared in units of m 2 takes the value 

2A(A + \h\) . 

Notice that the inflaton mass squared takes the analogous value for chaotic inflation in the h < — 1 regime changing 
D A while F{h) differs from G{h) given by eq. (|3.23p only by the sign of the 4 h 2 term. 

Recall that y ~ O(l) ~ h guarantees that g ~ O(10 -6 ) and A ~ O(10~ 12 ) without any fine tuning as stressed in 
ref. 0. 

In fig. [3l we plot 

w(x) as a function of ' 



8(h 2 + i) 2 8 VJ^TT 

for several values of h < 0. We see that the minimum of the potential 

Vvx+ =1 , \h\ 



V8Vh 2 + l Vh 2 + 1 

grows as \h\ grows. Similarly, the maximun of the potential at the origin 

yw(0) = F(h) 

8 {h 2 + l) 2 4 (h 2 + l) 2 

grows with \h\. 

New inflation is obtained by choosing the initial field x m the interval (0, x+)- The inflaton x slowly rolls down the 
slope of the potential from its initial value till the absolute minimum of the potential x+ ■ 

Computing the number of efolds from eq. (|2.8[) , we find the field x & t horizon crossing related to the coupling y and 
the asymmetry parameter h. We obtain by inserting eq. (|4.2p for w(x) into eq. (|2.8p and setting N[x] — N, 

y = z - 2 h 2 - 1 - 2 \h\ A + - \h\ (\h\ + A - </z) + 

3 



1 r 

— \h\ (A + \h\)A 2 log 



2F(/i)log[Vi(A-|/ l |)] , z^ V -x 2 - (4.5) 



z turns to be a monotonically decreasing function of y: z decreases from z = z + — (A + \h\) 2 till z — when y 
increases from y — till y = oa. When ^fz — > J~z+, y vanishes quadratically as, 



y = + 2 yi-y/zt) +0 {[yfr-y/*;] 



We obtain in analogous way from eqs. (|2.10p . (|2.1ip and (|4. 2|) the spectral index, its running, the ratio r and the 
amplitude of adiabatic perturbations, 

_ y z(z + 2feyi-l) 2 y 3z + 4feVi-l 

[F(h)-2z+lhzi/ 2 + z 2 ] 2 N F(h)-2z + §hz*/* + z* ' { ' } 
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FIG. 3: Trinomial New Inflation. The new inflation trinomial potential _ y ,^^.-> [eg. (I4.2M vs. the field variable 

different values of the asymmetry parameter /i = 0, —0.4, — 1, —20. We have normalized here the field variable and the potential 
by h-dependent factors in order to have a smooth \h\ S> 1 limit. 
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(4.7) 



(4.8) 



(4.9) 



A. The shallow limit (weak coupling) y — *• 

From eq. (|4.5p we see that in the shallow limit y — > 0, z tends to z+ = (A + |/i|) 2 . y(z) has its minimum y = at 
= z+. We find from eqs.(|4.5 ll -(|4.9 p . 



n s y =°l- 1-0.96 , r ^°A^o.l6, (4.10) 
A ~ nnnn. I A (S) ,2 V-0 ( M v ' 



^^-^,-0.0008 , lAn/^^^j A(A + , (4.11) 

which coincide with n s , f^ a k and r for the monomial quadratic potential. That is, the y — > limit is /i-independent 
except for |A^ d |. For fixed h and y — > the inflaton potential eq. (|4.2p becomes purely quadratic as we see from 
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eq. (|P|) : 

w( X ) V =° A(A + \h\) ( X - X+ f + 0{y/y) . (4.12) 

Notice that the amplitude of scalar adiabatic fluctuations cq. (|4.10j) turns out to be proportional to the square mass 
of the inflaton in this regime. We read this mass squared from eq. (|4.12[) : 2 A(A + \h\) in units of m 2 . The shift of 
the inflaton field by x+ has no observable consequences. For h = we recover the results of the monomial quadratic 
potential. 



B. The steep limit (strong coupling) y — > oo 

In the steep limit y — > oo, z tends to zero for new inflation. We find from eq. (|4.5 



y z =° -F(h) log z - q(h) -l + 0(Vz) , 



q(h) = 2 F(h) log (A - \h\) - - (h 2 + \h\ A) 8 A 2 log 



- 1-^ 



where 



q{h) is a monotonically increasing function of the asymmetry \h\ : < q(h) < oo for < \h\ < oo. 
Then, eqs.(|ii)] ) - (|4TF|) yield, 
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(4.13) 



(4.14) 



In the case ft, = we recover from the steep limit eas. (|4.13|l - (l4.14p the results for new inflation in the large y regime: 
we have -F(O) = 1 and q(0) = and eq. (|4.14p becomes, 



o 1 y 



a»x,h-»a 16 y _ x 
r = _ e 



y-»ljh->0 2 y 2 jfoj 

dlnfc TV 2 6 

|A (5) p,^ N 2 f M A 4 e^ +1 

1 feQrfl 12tt 2 VMf-J j/ 2 



(4.15) 



Notice that the slow- roll approximation is no longer valid when the coefficient of 1 /N becomes much larger than unity. 



C. The extremely asymmetric limit — > oo 
Eqs. (|4.5[ )- (|4.9p have a finite limit for \h\ — > oo with y and z scaling as /i 2 . Define, 

We have < Z < 4 for +oo > Y > 0. Then, we find for -> oo from eqs. (|4~5)) - ([49l) keeping Z and Y~ fixed, 

r = z-±VS-4-£ log* 16 



3 3 b 4 3 ' 

y z 2 (Vz-2f y zz-aVz 

n m-* z i + Z 2 ] 2 N m _ i zi + z 2 ' 1 ' ; 
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12 7T 2 \Mpi) Y 2 Z 2 (v r Z-2) 2 

In the — » oo limit, the inflaton potential takes the form 



w( x ) _ 32 l ,/F . y 4 



= lim -r^- = — --d—x l + — x 
|ft|-»oo h z AY 3 V 2 32 

This is a broken symmetric potential without quadratic term. Notice that the cubic coupling has dimension of a mass 
in eq. (|4.ip and hence this is not a massless potential contrary to the quartic monomial x 4 . 

In addition, eq. Q4.18p shows that for large \h\ one must keep the product \h\ M 2 fixed since this is determined 
by the amplitude of the adiabatic perturbations. We see from eq. Q4.18|) that in the \h\ — > oo limit, M = y/\h\ M 

becomes the energy scale of inflation. M ~ 10 16 GeV [eq. Q1.3|l ] according to the observed value of |A^ d |/iV displayed 
in eq. (!2.12|) . while M and m vanish as \h\ — > oo, 

M \h\^oo M 2 M 2 ihHoo 

M = — == = , ni = — — = -— — — — = . 

y^f M Pl \h\Mpi 



D. Regions of n 3 and r covered by New Inflation and by Chaotic Inflation. 

We see from eas. Q4.7p . (|4.10p and (I4.14p that new inflation for h < covers the narrow banana-shaped sector 
between the black lines in the (n s ,r) plane depicted in fig. [2] We have in this region: 

8 1 9236 
0<r<— , n s < 1 "' for oo > y > . (4.19) 

Chaotic inflation in the h < — 1 region covers the the even narrower complementary strip for -S. < r < t£ eq. (|3.26p . 
The zero coupling point y — > , r = ^ , n s = 1 — being the border between the two regimes. 

Chaotic inflation for — 1 < h < covers a wide region depicted by the red lines in the (n s ,r) plane [see fig. [2]. 
However, as shown by fig. [9l this wide area is only a small corner of the field z - asymmetry h plane. Since new 
inflation covers the banana-shaped region between the black curves, we see from this figure that the most probable 
values of r are definitely non-zero within trinomial new inflation. Precise lower bounds for r are derived from 
MCMC in eq.QU])- 



V. MONTE CARLO MARKOV CHAINS AND DATA ANALYSIS WITH THE TRINOMIAL 

INFLATION MODELS 

We performed Monte Carlo Markov Chains (MCMC) analysis of the commonly available CMB and LSS data using 
the CosmoMC program [ll|. For CMB we considered the three years WMAP data, which provide the dominating 
contribution, and also small scale data (ACBAR, CBI2, BOOMERANG03). For LSS we considered SDSS (DR4). 
We used the second release of WMAP likelihood code. In all our MCMC runs we did not marginalize over the SZ 
amplitude and did not include non-linear effects in the evolution of the matter spectrum. The relative corrections are 
in any case not significant [6j, especially in the present context. 

Given a theoretical model with several free parameters and a sufficiently rich set of empirical data, MCMC is a very 
efficient stochastic numerical method to approximately reconstruct the probability distribution for the actual values 
of those parameters. It is especially useful when the theoretical predictions are themselves stochastic in nature, as is 
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the case of the primordial fluctuations that seed the CMB anisotropies. CosmoMC is a publicly available FORTRAN 
program that performs MCMC analysis over the parameter space of the Standard Cosmological model and variations 
thereof. The spectral index n s of the adiabatic fluctuations and the ratio r of tensor to scalar fluctuations are two such 
parameters. In CosmoMC usage, one usually restrict the search to a subset of the full parameter space by imposing 
so-called hard constraints, whose number and type depend on previously acquired information both experimental and 
theoretical. 

Concerning statistical convergence tests, our MCMC data were collected in parallel Message Passing Interface (MPI) 
runs with 10 to 16 chains, with the 'R-l' stopping criterion (which looks at the fluctuations among parallel chains 
to decide when to stop the run) set to 0.03. Typically, this meant single chains with length of the order 10 5 , with 
the Gelman-Rubin convergence test [var(chain mean)/mean(chain var) in the second half of the chains] ranging from 
10~ 3 to 10~ 2 for the various MCMC parameters. 

We imposed as a hard constraint that r and n s are given by the analytic expressions at order 1 /N for the trinomial 
potential of the inflaton. Namely, r and n s are given by eqs. (|3.6p - (|3.8p for chaotic inflation and by eqs. (|4.6|) - (|4.8|) 
for new inflation. Our analysis differs in this crucial aspect from previous MCMC studies involving the WMAP3 
data [131 ]. As natural within inflation, we also included the inflationary consistency relation tit — — r/8 on the tensor 
spectral index. This constraint is in any case practically negligible. 

We allowed seven cosmological parameters to vary in our MCMC runs: the baryonic matter fraction ujb, the dark 
matter fraction uj c , the optical depth r, the ratio of the (approximate) sound horizon to the angular diameter distance 
9, the primordial superhorizon power in the curvature perturbation at 0.05 Mpc~\ A s , the scalar spectral index n s 
and the tensor-scalar ratio r. 

For comparison, we also run MCMC's within the standard ACDM model augmented by the tensor-scalar ratio r 
(ACDM+r model). That is, we treated n s and r as unconstrained Monte Carlo parameters, using standard priors. 
This analysis is indeed by now quite standard and good priors are available already in CosmoMC. 

We allowed the same seven parameters to vary in the MCMC runs for chaotic and new inflation. 

In the case of new inflation, since the characteristic banana-shaped allowed region in the (n s , r) plane is quite narrow 
and non-trivial, we used the two independent variables z and h in trinomial inflationary setup as MC parameters. 
That is, we used the analytic expressions at order l/N, eqs. (|4.6p and (|4.8p to express n s and r in terms of z and h. 
To be more precise, rather than z we used the appropriate normalized variable 

3i = l- — = 1 ■ - j (5.1) 

z+ (v^F+T+l/iQ 

z contains the field at horizon crossing and the coupling y. z+ stands for z at the absolute minimum of the potential. 
The variable z\ grows monotonically from to 1 as the coupling y grows from to oo [see eq. (|4.5I) and the lines right 
below it]. 

Concerning priors, we kept the same, standard ones, of the ACDM+r model for the first five parameters uj c , t, 9 
and A s ), while we considered all the possibilities for z\ and h, that is < Z\ < 1, < h < oo. 

In the case of chaotic inflation we kept n s and r as MC parameters, imposing as hard priors that they lay in the 
region described by chaotic inflation [see fig. This is technically convenient, since this region covers the major part 
of the probability support of n s and r in the ACDM+r model and the parametrization eqs. (|3.6[) - (|3.8|) in terms of the 
coupling parameters z and h becomes quite singular in the limit h — > — 1. This is indeed the limit which allows to 
cover the region of highest likelihood. 

The distributions for the field variable z and the asymmetry parameter h were obtained from the (n s , r) distributions 
by a numerical change of variables using eqs. ()3.6|) and (|3.8p and taking care of the corresponding Jacobian. We used 
a large collection of parallel chains with a total number of samples close to five million. In fact we did not need to 
explicitly compute the Jacobian. It was taken into account automatically since we used an uniform two dimensional 
grid in the z — h plane. Quite naturally in this approach for chaotic inflation, the most likely values of the cosmological 
parameters and the corresponding maximum of the likelihood coincide to those of the ACDM+r model. 

The priors on the other parameters where the same of the ACDM+r model and of new inflation. 

In all our MCMC runs we keep fixed the number of efolds N since horizon exit till the end of inflation. The reason is 
that the main physics that determines the value of N is not contained in the available data but involves the reheating 
era. Therefore, it is not reliable to fit N solely to the CMB and LSS data. 

The precise value of N is certainly near N = 50 [H, [13] ■ We take here the value N = 50 as a reference baseline 
value for numerical analysis, but from the explicit expressions obtained in the slow roll l/N expansion we see that 
both n s — 1 and r scale as l/N. 
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Therefore, decreasing or increasing N produces a scale trasformation in the (n r — l,r) plane, thus displacing the 
black and red curves in fig. [2]towards up and left or towards down and right. This produces however small quantitative 
changes in our bounds for r as well as in the most probable values for r and n s . We ran MCMC simulations with 
variable iV and imposing the trinomial new inflation potential. As a result, we found that N ~ 50 was the most 
probable value. A sharp fall in the likelihood happens for N < 50 and a significant decrease shows up for iV > 50. In 
summary, varying N has no dramatic effects and N ~ 50 turns to be the most probable value. The detailed analysis 
with variable N is beyond the scope of the present paper. 

We have not introduced the running of the spectral index dn s /d In k in our MCMC fits since the running [eq. (|2.Iip ] 
must be very small of the order 0(N~ 2 ) ~ 0.001 in slow-roll and for generic potentials Indeed we found that 
adding dn s /dlnk, as given by eq. (|3.7p or eq. (|4.7p . to the MCMC analysis yields insignificant changes on the fit of 
n s and r. 

On the contrary, when the running is introduced as a free par ameter, then the fit of n s and r gets worse and values 
for the running much larger than 0(N~ 2 ) ~ 0.001 follow pL Il3|. We think that the present data are not yet precise 
enough to allow a determination of dn s / dh\k. That is, adding further parameters to the fit (like the running) does 
not improve the fit and does not teach anything new. 



Our results for the new inflation trinomial potential eq. (|4.2p are summarized in figs. [4] and [6] 

In fig. 2] we plot the marginalized probability distributions (normalized to have maximum equal to one) of the most 
relevant cosmological parameters, which are the primary ones allowed to vary independently in the MCMC runs and 
some derived ones. The solid blue curves refer to the runs with the hard priors specific of trinomial new inflation. The 
dashed red curves are those of the ACDM+r model. As should have been expected from fig. [2j the really significant 
changes are restricted solely to n s and r. To provide further evidence of this, we compare in fig. [5]the joint probability 
(t, n s ) and (t, r) distributions of trinomial new inflation with those of the ACDM+r model. We recall that the optical 
depth parameter r is strongly correlated with n s . 

In the lower panels of fig. [6] we provide an enlarged version of the marginalized probability distributions for n s and 
r, together with their mean likelihoods. 

Notice the difference between probability distributions and mean likelihood distributions in the first four panels of 
both figs. [6]and[8] Probability and mean likelihood depend quite differently on the shapes and parametrizations of the 
allowed regions; if observational data were well concentrated within the allowed regions, both types of distributions 
would have a gaussian-like shape and the difference would be much smaller. With the current data, the joint probability 
distributions over the parameters of the trinomial potential are very far from gaussian, as can be appreciated from 



The theoretical constraints narrow the allowed region of parameters in a very nontrivial way. Otherwise, the 
parameters could cover a much wider region. Hence, the theoretically constrained distributions can hardly be and 
in fact are not gaussian. In the case of trinomial new inflation, the narrow banana-shaped region depicted in fig. [2] 
is responsible for the marked difference between marginalized probability distribution and mean likelihood for n s , as 
shown in fig. [SJ Besides, the banana-shape of trinomial new inflation in the n s — r plane, produces a spike in the left 
lower panel of fig. The sharp cut on the right is due to the theoretical upper bound on n s given by eq. (|4.I9p , the 
sharp rise on the left of the maximum is due to the marginalization over r. 

We find a lower bound on r for the trinomial new inflation model: 



A. 



MCMC results for new inflation. 



fig. m 



r > 0.016 (95% CL) 



r > 0.049 (68% CL) (new inflation) . 



(5.2) 



For n s we find: 



n s > 0.945 (95% CL) (new inflation) . 



The most likely values, which can be read off directly from MCMC data, are 



n s ~ 0.956 



r ~ 0.055 (new inflation) . 



The results for max(-log(likclihood)) are as follows 
ACDM+r model : 2714.038 



Trinomial chaotic inflation : 2714.038 



and 



Trinomial new inflation: 2714.153 
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FIG. 4: Comparison of the marginalized probability distributions (normalized to have maximum equal to one) of the most 
relevant cosmological parameters (both primary and derived) between Trinomial New Inflation (solid blue curves) and the 
ACDM+r' model (dashed red curves). 



At the point with most likely values for the cosmological parameters we see a small increase of 0.115 of 
— log(likelihood) for new inflation with respect to the ACDM+r model. This increase is only twice as much of 
the standard deviation of max(— log(likelihood)) over ten parallel chains, which is 0.062 in the case of the ACDM+r 
model. Therefore, this increase of the — log(likelihood) for new inflation is not significant. 

We recall that for trinomial new inflation there exist the theoretical upper limits: n s < 0.9615 . . . , r < 0.16 @ . 
That is, the most probable value of n s for trinomial new inflation is very close to its theoretical limiting value and that 
of r is not too far from it (see also fig. [6]). 

In the upper panels of fig. [5] we plot the marginalized probabilities and the mean likelihoods for the normalized 
coupling at horizon exit z\ and for y and h, that is the two parameters of the trinomial potential of new inflation. 
Their most likely values in our MCMC runs, 

Zi ~ 0.886, y ~ 2.01 and h ~ 0.266 (new inflation) , 

are not extremely significant in view of the shapes of probabilities and mean likelihoods in fig. [6] More appropriate 
quantifiers are the upper bounds 

y < 2.46 , h<1.2 with 68% CL , y < 6.35 , h < 4.92 with 95% CL (new inflation) (5.3) 

The z\ distribution decreases for small z\ (that is, where one approaches the quadratic monomial potential at y = 0) 
because such case yields a too large value for r, (r — > 0.16) with respect to the experimental data. For z — > 1 the z 
distribution decreases too because in this case one gets a too small spectral index n 8 . 

In conclusion, the most likely trinomial potentials for new inflation are almost symmetric and have moderate 
nonlinearity with the quartic coupling y of order 1. The \ — > — x symmetry is here broken since the absolute 
minimum of the potential is at \ ^ 0. 



22 






0.15 



FIG. 5: Upper panels: 95% and 68% percent contour plots of joint probability (r, n s ) distribution (left) and (r, r) distribution 
(right) in the ACDM+r model. Lower panels: the two joint distributions in trinomial new inflation. Recall that in this case 
n s < 0.9615 ... is a theoretical bound when N = 50. 

B. MCMC results for chaotic inflation. 

Our results for the trinomial potential of chaotic inflation, eq. (|4.2[) . are summarized in figs. \7\ [H] and [5] 
In fig. [7] we plot the marginalized probability distributions (normalized to have maximum equal to one) of the usual 
cosmological parameters as in fig.|4] Again, solid blue curves refer to the runs with the hard priors of chaotic inflation, 
while the dashed red curves are those of the ACDM+r model. As should have been expected from fig. [2J there are 
no really significant changes in any parameter. 

In the upper panels of fig. [5] we plot the marginalized probabilities and the mean likelihoods for the parameters of 
the trinomial potential of chaotic inflation, the normalized coupling at horizon exit z, the quartic coupling y and the 
asymmetry h. These are numerically calculated from those of n s and r (which are also reported in the lower panels 
of fig. |SJ), by means of eqs. (|3.6[) . (|3.8|) and (|3.24[) . The strong non-linearities of these relations are responsible for 
the peculiar shapes of the probabilities and mean likelihoods, and their marked relative differences, of the variables 
in the upper panels. 

Finally, in fig. [5]we depict the confidence levels at 12%, 27%, 45%, 68% and 95% in the (z, h) plane. The chaotic 
symmetric trinomial potential h = is almost certainly ruled out since h < —0.7 at 95 % confidence level. 

We see that the maximun probability is for strong asymmetry h < —0.9 and significant nonlinearity 0.5 < z < 1 
and 2.5 < y < 5. That is, all three terms in the trinomial potential w(x) do contribute to the same order. 

Notice that the range 0.5 < z < 1 corresponds for h — ► — 1 + to the region where the quartic coupling is quite 
significant: 4.207. . . < y < +oo according to eq. (|3.13p . 

The probability is maximum on a highly special and narrow corner in the (z, h) parameter space. This suggests: 

(i) the data force both the asymmetry as well as the coupling to be large. 

(ii) the fact that z ~ 1 implies large y means that we are in the nonlinear regime in z, that is strong coupling in \. 
This suggests that higher order terms may be added and will be relevant. 

(iii) If the preferred values are near the boundary of the parameter space, it could be that the true potential is 
beyond that boundary. The region of parameter space in which trinomial chaotic inflation yields r <C 1 is 
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FIG. 6: Trinomial New Inflation. Upper panels: probability distributions (solid blue curves) and mean likelihoods (dot-dashed 
red curves), all normalized to have maximum equal to one, for the values of the normalized coupling at horizon exit zi, of 
the quartic coupling y and of the modulus \h\ of the asymmetry of the potential. Lower panels: probabilities and mean 
likelihoods for the values of n s and r. Notice that here n s < 0.9615 . . . and r < 0.16 since these are the theoretical upper 
bounds for n s and r in trinomial new inflation with N = 50. Trinomial new inflation exhibits a lower bound for the ratio: 
r > 0.016 (95% CL), r > 0.049 (68% CL). 



very narrow and highly non-generic and corresponds to an infiaton potential contrary to the Landau-Ginsburg 
spirit since adding higher degree terms in the field to the trinomial potential can produce relevant changes. This 
means that the fit for chaotic inflation can be unstable while for new inflation the maximun probability happens 
for a moderate nonlinearity and therefore will not be much affected by higher degree monomials. Finally, the 
r — > regime is obtained near the singular point z = 1, h = — 1 where the infiaton potential vanishes. 

(iv) The MCMC runs appear to go towards the h = — 1 limiting chaotic potential exhibiting an inflexion point [sec 
fig. [T] which is at the boundary of the space of parameters. This may indicate that the true potential is not 
within the class of chaotic potentials. Since the runs go towards the maximal value of the asymmetry parameter 
h where the infiaton potential exhibits an inflexion point, this supports the idea that the true potential must 
definitely break the \ ~~ * ~ X symmetry, as new infiaton potentials do spontaneously. This favours again new 
inflation since the best fit to the trinomial new inflation potential corresponds to small or zero asymmetry 
parameter h. The spontaneous breakdown of the symmetry seems enough to obtain the best fit to the data 
without the presence of an explicit symmetry breaking term | x 3 [ see eq. (|4.2p ]. 



C. Conclusions. 



The MCMC results for chaotic and new inflation presented above show that: 

(i) symmetry breaking is begged by the data. Namely, in chaotic inflation the data ask for a strongly asymmetric 
potential since we find the maximun probability when h is near the extreme value h = — 1. In chaotic inflation, 
the h = case is almost certainly ruled out as noticed above. In new inflation the symmetry is spontaneously 
broken by construction since the absolute minima is at \ ^ 0. Hence, in one way or the other the data request 
a breakdown either explicit or spontaneous of the \ ^ ~X symmetry. 

New inflation naturally satisfies this requirement since the spontaneous symmetry breaking can alone reproduce 
the data with a moderate nonlinearity. 

(ii) Adding higher order terms to the trinomial new inflation potential will not affect the trinomial fit since it 
corresponds to a moderate nonlinearity z\ = 0.886. On the contrary, adding higher order terms in trinomial 
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FIG. 7: Comparison of the marginalized probability distributions (normalized to have maximum equal to one) of the most 
relevant cosmological parameters (both primary and derived) between Trinomial Chaotic Inflation (solid blue curves) and the 
ACDM+r model (dashed red curves) . 



chaotic inflation can change the trinomial fit since we find a significative nonlinearity 0.5 < z < 1 with the 
trinomial chaotic potential. Therefore, in chaotic inflation, each time a new term of higher degree is added to 
the potential, it changes significatively all the coefficients in the potential. This is against the Landau-Ginsburg 
spirit which is the framework of the effective theory of inflation.] 

(iii) fnflaton potentials of degree higher than four have been recently investigated in ref. [ll|. It is shown in ref. 
[l5l | that binomial potentials with quadratic and \ 2 " terms reproduce the data in domains which are larger for 
n = 2 than for n > 2. In addition, for new inflation it is shown that the crucial region to reproduce the data 
is the neighbourhood of the absolute minima of the potential. That is, only the quadratic and somehow the 
quartic terms around x+ are relevant to fit the data. 

In summary, the MCMC analysis of the CMB+LSS data indicates that new trinomial inflation is the best model. 

Trinomial new inflation gives a very good fit to the data (CMB+LSS) with a moderate nonlinearity and small 
or zero asymmetry parameter. That is, the potential which best fits the data is the spontaneously broken 
symmetry potential eq. (|4.2p with h = and y ~ 2.01 . . . 

, 2 



W (X) = ( X 2 



I) 



Moreover, we have for this potential, when N = 50, 

n s ~ 0.956 , r ~ 0.055 (new inflation) . 

(iv) More generally, we have a lower bound on the ratio r within the trinomial new inflationary potentials given 
by 

r > 0.016 (95% CL), r > 0.049 (68% CL) (new inflation) . 
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